Competing Order in the Mixed State of High Temperature Superconductors 
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We examine the low temperature behavior of the mixed state of a layered superconductor in 
the vicinity of a quantum critical point separating a pure superconducting phase from a phase in 
which a competing order coexists with superconductivity. At zero temperature, we find that there 
is an avoided critical point in the sense that the phase boundary in the limit B — > does not 
connect to the B — critical point. Consequently, there exists a quasi- ID regime of the phase 
diagram, in which the competing order is largely confined to ID "halos" about each vortex core, 
and in which interactions between neighboring vortices, although relevant at low temperature, are 
relatively weak. 

PACS numbers: 



Whereas in many well understood metallic compounds 
over a broad range of compositions and temperatures, 
the only two phases encountered are the normal (Fermi 
liquid) and superconducting phases, in the cuprate high 
temperature superconductors, and other highly corre- 
lated electronic systems, there are many ordered phases 
which appear to compete and sometimes coexist. In 
addition to the uniform d-wave superconducting state, 
compelling evidence exists of ordered antiferromangetic 
(Neel), unidirectional charge density wave ("charge- 
stripe"), unidirectional, colinear, incommensurate spin- 
density wave ( "spin-stripe" ) phases and those with coex- 
isting superconducting and stripe order jjj]. Preliminary 
evidence also exists of possible d-density-wave or stag- 
gered flux order ||, ||, electron nematic Q, d + id or 
d + is superconducting order |^|, ^j, and various other 
phases in which more than one of these orders coexist. 

One rather direct way to look for competing order pa- 
rameters in a system which is globally superconducting 
was recently proposed by Zhang and coworkers [Q, ||] and 
others |], [n| [ll]]. The idea is that if the superconducting 
order is only slightly favored over a competing order, then 
where the superconducting order is suppressed in the core 
of a vortex, the competing order will be manifest. Indeed, 
recent neutron scattering experiments of Lake et. al. 
on La2- :r Sr 2 ,Cu04 have revealed a strong enhancement 
of "spin-stripe" order at low energies produced by mod- 
est magnetic fields. Similar results have been obtained in 
La2Cu04 by Birgeneau et. al. |l5|] . Moreover, in scanning 
tunnelling microscopy (STM) studies by J. Hoffman et. 
al. of near optimally doped P^S^CaC^Os+a in a 
7T magnetic field, large induced "halos" about each vor- 
tex core have been imaged where the density of states 
is modulated with a spatial period (4a) equal to that 
expected for "charge-stripe" order. Additional evi- 
dence that there is a substantial degree of local charge 
stripe order in Bi2Sr2CaCu20s+<5 comes from the similar 
patterns of density of states modulation observed in the 



zero-field STM studies of C. Howald et. al. jL6j on the 
same material. (This interpretation is, however, being 
challenged in a forthcoming paper by Hoffman et q/p7|.) 

The notion of a competing order developing an expec- 
tation value in a vortex core is basically a mean-field no- 
tion. However, in all the examples mentioned above, the 
competing order is associated with spontaneously broken 
symmetry, so that fluctuation effects may fundamentally 
alter the physics. Specifically, in a planar system, the 
vortex core is a finite size system, and so cannot support 
a spontaneously broken symmetry, while in three dimen- 
sional superconductors the vortex core is a one dimen- 
sional system, so cannot exhibit any symmetry breaking 
except at T — 0, and there only for discrete symmetries. 



I. PHASE DIAGRAM 

Our principal results are summarized in the schematic 
zero temperature phase diagram (Fig. |l|) for a layered 
system in which superconducting order, with an order 
parameter denoted by ^ , competes with another type of 
order, whose order parameter is denoted by <p. The two 
axes represent the magnetic induction, B, and a control 
parameter, "a," such as pressure or doping concentra- 
tion, with the convention that increasing a disfavors <f> 
order. It is assumed that at B = there exists a con- 
tinuous quantum phase transition at a — a c (the heavy 
circle) separating a pure superconducting phase from a 
phase with coexisting (j) and superconducting order. 

While the considerations here are rather general, it is 
useful to have specific realizations in mind. For us, the 
most important case is that in which <f> represents |Is[ ] 
stripe orientational ("nematic") order. In the absence of 
crystal field effects, the orientation of the stripes is arbi- 
trary and can be parametrized by an angle < < n. 
Hence, in this case cb would then be a director, a head- 
less vector, which in two dimensions can be represented 
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FIG. 1: Schematic zero temperature phase diagram of a lay- 
ered superconductor as a function of a control parameter a 
(described in the text) and the magnetic induction B, i.e. 
the actual magnetic field which penetrates the system. Here 
*P and 4> denote, respectively, the expectation values of the 
superconducting and competing order parameters and it is 
assumed that for B = there is a continuous transition at 
a — a c between a pure superconducting phase and a phase 
with a coexisting and cf> order. The solid lines denote phase 
boundaries and the other lines are crossovers (all are described 
in the text). In particular, a c (B) represents the boundary be- 
tween the pure superconducting and the phase with coexisting 
<fi and superconducting order, and a± = lim_g^o at c (B) marks 
the point at which the <f) "halo" about an isolated vortex-line 
undergoes a transition from a quantum disordered state (for 
a > ai) to an ordered or quasi-ordered state for a < a.\. 
(af is the mean-field value of et\.) 



by a complex scalar field defined so that <fi = \(f>\ e z2e 
corresponds to stripes lying along a preferred direction 
eo = xcos9 + ysm9. (Note that 9 = and 9 = tt are 
physically equivalent.) In this limit, <p nas a continuous 
XY symmetry. However, in a crystal with appropriate 
point-group symmetry |l9[ there are only two preferred 
stripe orientations. In this case, </> reduces to a real scalar 
field reflecting the Ising character of the symmetry break- 
ing. 

If, on the other hand the relevant stripe order is mag- 
netic (incommensurate SDW) order, in addition to ori- 
cntational symmetry the relevant broken symmetry is 
spin rotational invariance, so for each stripe orientation, 
</> is a three component real vector field corresponding 
to Heisenberg symmetry. We shall see that this case 
is slightly different than the lower symmetry situations. 
Stripe states can also spontaneously break translational 
symmetry, but since a vortex core explicitly breaks this 
symmetry in any case, issues of spatial symmetry break- 
ing are more subtle, and will be discussed elsewhere [ p0[ . 

In the phase diagram shown in Fig. [j], the solid lines 
represent actual phase boundaries, and the broken lines 
are crossovers. The most striking feature of this phase 
diagram is the presence pl[ of an "avoided critical point," 



i. e. the phase boundary a c (B) has a discontinuity at 
B = 0: 

lim a c (B) ^ a c (1) 

This discontinuity in the phase boundary is a conse- 
quence of the fact that the magnetic field is a singular 
perturbation. 

The assumed competition (mutual suppression) be- 
tween superconducting and <p order (the latter assumed 
to be essentially uncoupled to B) is seen in the fact that 
the critical line, a = a c (B) separating the pure and co- 
existence phases is an increasing function of B up to 
the point at which superconductivity is completely sup- 
pressed by a magnetic field in excess of B = H C 2- The 
dashed line marks a crossover from quasi one dimensional 
to fully three dimensional order (where <j) is more or less 
uniform in space). It essentially coincides with the phase 
boundary derived by Demler et. al. p~T[ ] for the case in 
which interactions between layers are neglected^]. The 
dotted line represents a crossover associated with the 
mean-field phase boundary; just to the left of this line, 
there is a "halo" of the competing order surrounding each 
vortex core, but each halo fluctuates essentially indepen- 
dently, and there is no true <ft order. 

The arguments that lead to this phase diagram consti- 
tute the bulk of the present paper. The reader should be 
warned that there are some subtleties (which we will dis- 
cuss below) associated with one or another specific type 
of competing order that can affect the shape, and even 
the topology of the phase diagram, as does even weak 
quenched disorder. For simplicity of discussion, for most 
of the paper we will assume the extreme type II limit, in 
which the London penetration depth, A = oo, although 
this assumption is not necessary. Some of the experimen- 
tal consequences of this phase diagram are discussed in 
the final section of this paper. 

II. THE BASIC PHYSICS 

In this section, we sketch the basic physics that leads 
to the phase diagram in Fig. [I]. All actual derivations 
are deferred to later sections. For simplicity, let us first 
consider the case in which the broken symmetry associ- 
ated with the competing order is Ising-like, i. e. is a 
real scalar, and there is a symmetry under <fi — > —cj). 

A. The Ising Case 

We start our discussion by considering the structure of 
a single, isolated vortex in the uniform superconducting 
state j23|, i. e. a > a c . If a is large, the structure of the 
vortex is unaffected by the proximity of the (f> ordered 
phase, but if a is sufficiently close to a c , then at mean- 
field level, the suppression of \& in a vortex core region 
of size £o (the superconducting coherence length) will re- 
sult in a halo with radius equal to the critical correlation 
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length, £0 3> £0) i n which <p 7^ 0. An important aspect 
of the structure of the vortex, first emphasized in this 
context by Demler et. al. [0, is that the magnitude of 
the superconducting order parameter does not return to 
its bulk value exponentially, but rather, for a vortex at 
the origin and for r>(o, has a power-law form 



tO X<f>(T c )Jinte 



1. This leads to the estimate 



l^l 2 -^ [l-(Wrf + 



(2) 



This 1 jr 2 fall-off is a necessary consequence of the slow 
decay of the superfluid flow around an isolated vortex, 
and results in a somewhat larger halo of 4> order around 
a vortex than would otherwise occur. In any case, £^ ~ 
\a — a c \~ u diverges with a critical exponent i/asa-> a c . 
Since v is the quantum critical exponent of a system in 
d = 3 spatial dimensions, it presumably takes its mean- 
field value, v = 1/2 (up to logarithmic corrections-to- 
scaling). 

Within a single superconducting plane the halo is a 
finite size system, and so cannot actually support a bro- 
ken symmetry state, i. e. quantum fluctuations will cause 
the system to tunnel between the positive and negative 
4> states with a matrix element h = h(a). Thus, in the 
absence of inter-plane couplings, there can be no true 
<f> order established until B is large enough, i. e. the 
vortex density is large enough, that there is significant 
coupling between neighboring vortices. Because the size 
of the halo increases with increasing we expect that 
h vanishes in the limit £^/£o 00 ■ We derive an explicit 
expression for this, below in Eq. |2^, from the Landau- 
Ginzburg theory of competing orders. However, even if 
the microscopic coupling between the order parameter <j> 
on neighboring planes, Jo, is weak, the large size of the 
halo implies a large effective coupling, J, between the 
halos on neighboring planes, with J ~ Jo^ 2 ^ 2 ,! where 

(j) 2 is the mean squared value of <f> in the vortex halo. 
From the Landau- Ginzburg theory, this coupling is seen 
to diverge as £,<p/£,o — ► 00, although only logarithmically, 
J oc J log(^/^o)j as shown in Eqs. ([l3]) and (|l9|). 

Thus, along an isolated vortex, the low energy fluctu- 
ations are equivalent to an effective transverse-field Ising 
model with Ising coupling J and transverse field h. This 
model has an ordered ground-state provided J/h > 1. 
Since as a approaches a c , the limit J/h — > 00 is real- 
ized, it follows that for a near a c , an isolated vortex line 
has an ordered ground-state! Moreover, in this limit, 
no matter how small B, the inter- vortex coupling can 
never be ignored, and indeed leads to a finite transition 
temperature. We can estimate this transition temper- 
ature using standard methods |?4j of quasi-one dimen- 
sional systems: We estimate the inter-vortex coupling to 

be Jinter ~ exp(-r/^) — exp(-A/ \J Bt, 2 ), where r is 

the spacing between vortices and A is a number of order 
1. At low temperatures, we estimate the susceptibility of 
an isolated vortex line to be that of the ID Ising chain, 
X<t>(T) ~ exp(A'J/T). Finally, we estimate T c according 



T C (B) ~ fyjVB 



(3) 



£?V2| a _ ac |-V2 log [ ac /| a _ ac |] 



where in the second expression we have adopted the 
Landau-Ginzburg estimate of J. 

Reversing the present logic, it is clear that with increas- 
ing a, we will eventually encounter a condition in which 
h > J. Here, the ground-state of the isolated vortex line 
is quantum disordered, and hence even at T = 0, there 
is no <fi ordering until a critical field strength, B c (a), is 
exceeded, such that the interactions between neighbor- 
ing vortices is strong enough to induce ordering. The 
critical value a c (B) — > ol\ in the limit B — > marks 
the quantum critical point of the isolated vortex, where 
J(a>i)/h(a>i) = 1. (For small \a— ai\, the phase boundary 
is non-analytic |2a] - a subtlety which we have neglected 
in sketching Fig. hi) 

The remaining phase boundaries that occur at larger 
values of B are determined by the obvious and conven- 
tional physics of competing orders. We now discuss the 
physics of the two crossover lines shown in Fig. |]. 

Where the <j> halos about the vortices start to over- 
lap strongly, there is a crossover from a quasi-one dimen- 
sional regime in which the magnitude of <fi is substantially 
inhomogeneous, to a regime where the variations of (f> are 
relatively small. (Indicated by the dashed line in Fig. |I|.) 
At first guess might be that this crossover occurs when 
the spacing between vortices is of order £^ or in other 
words when B* ~ \a — a c \ 2v = \a — a c \ (up to logarith- 
mic corrections to scaling). In fact, as shown by Demler 
et. al. j 1 1 , this crossover occurs at somewhat smaller B 



due to the slow recovery of \& away from the vortex core. 
If we consider the case in which <jj is homogeneous, then 
it is as if there were a magnetic field dependent reduc- 
tion of the effective a G ff = a — */\B\ log \H C 2/B\ where, 
as we will see in the next section, 7 is proportional to 
the coupling strength between the two order parameters. 
Consequently B* ~ \a — a c \/ log \a c /(a — a c )\. (Pre- 
cisely at a = a c , the same line of reasoning leads to the 
conclusion that T c ~ \Blog(H C 2/B)\" z , where z is the 
dynamical scaling exponent, and hence depends on the 
dynamics of <p.) 

For large enough a > , the competing order is suf- 
ficiently disfavored that, even at mean field level, no halo 
of non-zero (j) is induced about the vortex core. Thus, 
for small B and a c (B) < a < af \ there are substantial 
local <fi correlations along each vortex, but the quantum 
fluctuations are sufficiently large that the coupling be- 
tween vortices can be neglected, and no true long-range 
<j) order develops. For a > a^, there are no substan- 
tial 4> fluctuations induced by the presence of vortices. 
(The value of a = o^ ' at which a non-zero value of <f> 
appears at mean-field level is computed from Landau- 
Ginzburg theory in Eq. (|l^); roughly it is the point at 
which £0 — £o-) 
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We note, in passing, that even in 2D the Ising case is 
quite different than than Hciscnberg case considered by 
Demler et aZfff]]. In particular, as discussed in footnote 
pU, in the Ising case the transition to the coexistence 
phase occurs when the spacing between vortices is para- 
metrically large compared to Consequently, although 
linis^o ct c (B) — a c , for a close to a c (B) and B small, <p 
is not spatially uniform but rather is strongly peaked in 
halos about individual vortex cores. 



B. XY and Heisenberg cases 

If <j) is not an Ising variable, but has higher symme- 
try, the above considerations are somewhat modified. In 
the case of an XY variable, the physics of the isolated 
vortex is equivalent to the well-known physics of the ID 
quantum rotor model. Again, there is a single site {i.e. 
intra-plane) term, which can be characterized by the en- 
ergy gap, h, between the ground-state and first excited 
state of a single rotor; manifestly, h is proportional to 
the inverse moment of inertia of the rotor and large h 
favors a quantum disordered state. Ordering, however, 
is again promoted by an "exchange" interaction, J, be- 
tween neighboring "sites," i. e. neighboring planes. As 
in the Ising case, we estimate the dependence of J on a 
from the Landau-Ginzburg treatment, below, and with 
the same result. For a continuous symmetry, h does not 
involve tunneling, and so its dependence on a is much 
weaker than in the Ising case; indeed, we will see from 
the Landau-Ginzburg treatment that h has the inverse 
dependence as J, i. e. h(a) cx 1/J(a). 

There is still a quantum disordered phase possible if 
h/J is sufficiently large. However, at smaller h/J the 
ordered phase of the Ising chain is replaced by a confor- 
mally invariant (power-law) phase in the XY chain. The 
susceptibility in this phase is a power-law in tempera- 
ture, X<f>(T) ~ J^ 1 (J/T) K , where K{J/h) is an increas- 
ing function of J/h. So long as K > 0, this susceptibility 
still diverges as T — > 0, so the topology of the phase di- 
agram is similar to that in the Ising case. However, now 
a.\ is determined implicitly from the relation K(a±) = 0. 
For ct\ > a > a c , where K(a) > 0, we can use the same 
inter-vortex mean-field theory to estimate the ordering 
temperature: 



T c ~ Jexp[-A/(K£+VB)] 



(4) 



i. e. T c rapidly becomes immeasurably small at small B. 

The classical Heisenberg ferromagnet has an ordered 
ground state, even in one dimension. However, the 
0(3) quantum rotor model, which represents the low 
energy theory of the quantum Heisenberg antiferromag- 
nets, does not. In the absence of Berry-phase terms, 
the one-dimensional 0(3) chain possesses [^6| a Hal- 
dane gap, and hence a non-divergent zero temperature 
susceptibility. This means that the topology of the 
phase diagram is different for the Heisenberg case, and 
that there is no avoided critical point. However, since 



the Haldane gap vanishes exponentially for large J/h, 
Ejjaidane ~ exp[— A" J / h] , it follows from simple scaling 
arguments that x<j>( T = °) ~ exp[2 A" J/h]. Thus, fol- 
lowing a line of argument similar to the one which in the 
Ising case led to Eq. ([IJ), inter- vortex mean field theory 
leads to an estimate for the critical B, 



B c (a) ~ [&J/h]- 



|(«A* C )-1| 
log 2 |(a/a c )-l| 



(5) 



where in the second expression we have adopted the 
Landau-Ginzburg estimates of h and J. 

It is also possible to imagine that a single vortex cor- 
responds to a half-integer spin chain, in which case there 
is a Berry's phase and consequently no Haldane gap. In 
this case, the situation is essentially equivalent to the XY 
case, with the susceptibility exponent K — 1 (up to log- 
arithmic corrections). However, it seems to us that since 
the effective Heisenberg model in the present case is not 
sharply defined on the lattice scale, fluctuation effects 
are likely to smear out the subtle interference phenom- 
ena responsible for the special behavior of half-integer 
spins, leaving us with the physics of the rotor described 
above. 



C. Further subtleties 

There are still other subtleties to worry about. In this 
entire discussion we have assumed that the vortex texture 
of ty(r) does not lift the symmetries which are sponta- 
neously broken by the ordering of cf>. However, where one 
of those broken symmetries is translation invariance, the 
presence of a vortex core is an explicit symmetry breaking 
field. If this effect is significant (as it may well be in the 
case of stripe order), it greatly complicates the analysis. 
A related issue is that we have assumed that there is no 
frustration of global <fi order which arises from the form of 
the vortex lattice and the nature of the coupling between 
neighboring halos. We believe that, in the absence of the 
just mentioned symmetry breaking terms, this assump- 
tion is reasonable, but it is manifestly unreasonable in 
their presence. Finally, especially in the regime where the 
physics is quasi one dimensional, all results are likely to 
be extremely sensitive to even tiny amounts of quenched 
disorder. We have not fully explored the implications of 
any of these further problems. 



III. LANDAU-GINZBURG THEORY 

While most features of the problem are largely deter- 
mined by considerations of order parameter symmetry, 
it is pedagogically useful to make them explicit by con- 
sidering the Landau-Ginzburg treatment of two compet- 
ing order parameters. We are interested in ground-state 
properties, so we must ultimately analyze a D+l dimen- 
sional quantum action, where D=3 is the spatial dimen- 
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sion. However, we will be analyzing this action semi- 
classically, in the sense that we will first consider time- 
independent field configurations which minimize the ac- 
tion, and then analyze quantum fluctuations about this 
classical ground-state. Thus, we start by considering only 
the classical (static) Landau-Ginzburg free energy den- 
sity functional in a single plane of a layered system 



2 

2 



(6) 



where 'J, a complex scalar field, is the superconducting 
order parameter and 0, which may have multiple com- 
ponents, represents the competing order parameter. In 
the present paper, we will only focus on "competing" or- 
ders in the sense that we will always assume that 7 > 0. 
Indeed, we assume that 7 is not small - if 7 is small it 
means that the two order parameters hardly interact, as 
happens, for instance, in conventional superconductors 
when and \& originate from different pieces of the Fermi 
surface. The free energy must be invariant under a global 
U(l) transformation, fy(r) — > e lc * 0> l'(r) due to gauge in- 
variance, and, depending on the nature of the competing 
order, under an additional set of global transformations, 
<j){r) — ► g4>(f) where g are elements of an appropriate 
coset space. For instance, if 4> corresponds to Neel order, 
then g € SU(2)/U(1). In Eq. |[ the superconducting 
contribution to T is of the usual form, 



2 K i 



2e 



A)^\ 2 



\J>2 



(7) 



where A is the vector potential, £0 = y/ Ko /2^q is the 
coherence length, and here, and elsewhere, . . . refers to 
higher order terms in powers of the order parameters. 
The competing order is governed by 



K<t> 



2 |V0| 2 + ||0| 



a, ,,2 , 1, 



(8) 



The mean-field solution is obtained by solving 
the Landau-Ginzburg (LG) equations, 8J-/5^\q, = 
8J- /5<j)\^ = 0. For the most part, we will focus on states 
deep in the superconducting phase, where we can treat 
*&(f) as a given function, leaving us with the task of com- 
puting <p. By symmetry, <fi — is always a solution; where 
a non-trivial solution exists, it leads to a condensation 
energy: 



E, 



cond 



= / dr[F(*,t)-F(*,0)] = - 



In the spatially uniform case, \&(f) 
a non-trivial solution 

a < a r — 



(9) 



^0) there is 



r) = yj{a c - a) /{I - 7 2 ) for 
7 'J 2 ,, while <f> = for a > a c . In the cf> 
disordered phase, it is easy to see that any inhomoge- 
neous solution of the LG equations will decay exponen- 
tially with a correlation length ^ = y/n^Ka — a c ). 

In the case of a single vortex at the origin, we look 
for a solution of the form <P(r) = e lv |^(r)| where <p is 



the azimuthal angle. While the exact form of ^(r)! is 
somewhat complicated, at large r it is easily seen to be 
|¥(r)| = * [1 - (l/2)(f /r) 2 + O((£o/r) 4 )]. (Of course, 
if magnetic screening is taken into account, |\t(r)| ap- 
proaches ^0 exponentially for r > A, with the London 
penetration length; in the high temperature supercon- 
ductors, A is very large, so it is reasonable to approxi- 
mate it as infinite.) comment at various stages on the 
effect of finite A on our results.) Since none of our results 
depend critically on the short distance (r < £0) behavior 
of this solution, we will adopt the approximation 



|*(r)| 2 = - (£o/r) 2 ] for r > £0 



(10) 



and ^(r)! 2 = for r < £q- Thus, the LG equation for 
an isolated vortex is 



dx 2 



ld_ 

x dx 



£0 



(11) 



where x = r/£ , 7 = 2^ 2 ^ 2 /n l f > and 4>(r) = 
(y/K$/£o) <t>{%)- This equation is valid for x > 1; for 
x < 1, the term 7/x 2 is replaced by 7. Manifestly, at 



large x, <j>(x) 



-1/2 



exp(— x(£o/^)) falls with the ap- 



propriate Ornstein-Zernicke form, so where a non-trivial 
solution exists, the condensation energy is always finite. 
(As expected, in unsealed units, the vortex halo has a 
radius £<£.) 

At the critical point, a — a c (i.e. £o/£<£ = 0)) it is easy 
to see that the solution of Eq. (O) is 



5 ) = yl + 7 x 1 for x > 1. 



(12) 



(In the screened case, for x > (A/£o), the solution is of the 
same form, but with VI + 7 — ► 1.) It thus follows that 
for large £0/^0 > 4> looks critical for a large intermediate 
range, £^ > r > £o- From this, it follows that, asa-> a+, 



dr\0\ 2 
df\0\ 4 



27r^(l + 7)log|AC /£ o | 



(13) 



where A and A' are numbers of order 1 . The first of these 
results is the essential ingredient leading to the existence 
of an avoided critical point in this problem. Note that 
the long-range tails of the vortex profile contribute to, 
but are not essential to this result; even if the vortex is 
screened, (f> at the critical coupling still has a 1/r form, so 
the second moment still diverges logarithmically as the 
critical coupling is approached, although with a some- 
what different pre-factor. 

For a not so close to a c , we expect a critical value 
of a — such that for a > there are no, non- 
trivial solutions to the LG equations in the presence of an 
isolated vortex. It is relatively straightforward to prove 
that a necessary and sufficient condition for the existence 



6 



of a non-trivial solution is that the quadratic kernel in T 
have at least one negative eigenvalue, i. e. that 



dx 1 



ld_ 

x dx 



■ 7 |*(r)P 



Mr) = e<MO (14) 



has a solution with e < 0. Thus, a\ is defined implicitly 
from the condition e = 0. For 7=1 and the approximate 
form of the vortex profile introduced above, 



„(°) 



^=-0.75-^ or & = 1.97£o. 
So 



(15) 



Similar estimates can be obtained for any positive 7. 

Finally, it is clear from the asymptotic form of the 
single vortex solutions that for dilute vortices, the effects 
of coupling between vortices are of order 



Jinter ~ exp[-i?/^] 

where R 3> £^ is the spacing between vortices. 

IV. THE EFFECTIVE ACTION 



(16) 



Our goal is to obtain explicit expressions for the cou- 
plings in the effective Hamiltonian, discussed in Sec. II, 
above. In particular, wherever there is a non-trivial solu- 
tion to the mean-field equations, there are clearly a family 
of equivalent solutions, and important fluctuations which 
could potentially restore the broken symmetry are those 
which carry the system from one solution to another. 

We begin with the case in which the order parameter 
(f> has Ising symmetry. In the Ising case, the classical 
ground-states are <j> = where a — ±1. Thus, at 
low low temperatures, the relevant states are of the form 
4>j{f) = <jj4>(r). The effective Hamiltonian has matrix 
elements, denoted below by h, connecting these classi- 
cal states with each other via tunneling processes ( "spin 
flip" ) . However, while a detailed calculation of these ma- 
trix elements requires additional analysis to determine 
the tunneling paths that locally permit the system to 
fluctuate from one ground-state to the other, the end re- 
sult is quite simple: the resulting effective Hamiltonian 
must be of the form of a transverse-field Ising model 



H. 



off 



(18) 



To complete the quantum description of the order pa- 
rameter fluctuations in the presence of an isolated vortex- 
line, we define an effective Euclidean action which in- 
cludes inter-plane couplings and the simplest possible 
quantum dynamics. We will consider the dynamics only 
of the slow modes of the order parameter <fi. However we 
will ignore the effects of the (possibly interesting) cou- 
pling to superconducting quasiparticles and other low 
energy degrees of freedom - at least in a BCS d-wave 
superconductor they have a vanishingly small density of 
states and their effects on the dynamics of the order pa- 
rameters of interest here is likely small. In addition, since 
the phase mode of the superconducting order parameter 
* decouples at long- wave-length from the rest of the de- 
grees of freedom, it is reasonable to suppose that quan- 
tum fluctuations of the superconducting order parameter 
can be integrated out to produce a small renormalization 
of the effective parameters; we thus consider only static 
configurations of 'J', and omit any dependence on . For 
all the cases of competing order <f> that we are discussing 
here, the effective Euclidean action takes the form 

S[^}= Y^J^dr J dr[{M/2M 2 (17) 

-J (4^- +1 +h. c.)} 

where r is a 2D vector in each plane, j labels the planes, 
4>j{r) is the order parameter field in plane j, and we only 
consider the case in which = W(f) is independent 

of layer index, i.e. when we consider a vortex core, we 
assume the vortex line is static and precisely perpendic- 
ular to the layers. Other forms of the <f> dynamics can be 
considered, but we believe that, for the most part, the 
results will not be qualitatively different. 



where J is given by 



J/2J 



df\</)(r)\ 2 



(19) 



In Eq. [18|, h is the tunneling matrix element which there- 
fore depends exponentially on the product of the effective 
mass M* times the barrier height -E C ond- Since the effec- 
tive mass is renormalized by precisely the same factor, 
M*/M = J/2 J , it follows that 



logh~-A"y/M*\Ec 



log[W£o] 



(20) 



where the prefactor A" is a constant determined by de- 
tails of the tunneling process. This implies that this 
problem approaches the classical ID Ising- ferromagnet 
as a — > a c . 

This sort of analysis becomes simpler in the case in 
which a continuous symmetry is broken by (j). For in- 
stance, consider the XY and Heisenberg cases in which cf> 
is, respectively, a two or three component vector, while <j), 
which is the solution of the LG equations, can be taken 
to be a scalar. Then, the full set of mean-field solutions 
can be written as 4>{r) = Cl(f)(r), where fi is a unit vec- 
tor in order-parameter space. Now, for the case of weak 
inter-layer coupling and large "mass," M* , we can ignore 
"amplitude" fluctuations, and focus our attention on the 
soft Goldstone modes by restricting our attention to field 
configurations of the form 



so that the effective Hamiltonian is 



(21) 



(22) 
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where Ly is the angular-momentum conjugate to fi, and 
J is given by the same expression, Eq. ( ]19| ) , as in the Ising 
case. This is the Hamiltonian of a ID array of quantum 
rotors, as discussed previously. Clearly, it follows from 
Eq. Jl3| ) that asa^a c , both J and M* diverge, making 
the system more and more nearly a classical ferromagnet. 

These results flesh out the general physical arguments 
made in the beginning of this paper. The effect of am- 
plitude fluctuations have not been included in any of the 
present considerations - since ultimately we are dealing 
with a 3D quantum system, we believe they are gener- 
ally less important than the "phase" fluctuations we have 
explicitly treated. We have not analyzed the further fluc- 
tuation corrections to test this supposition. 

V. CONCLUDING REMARKS 

When there is more than one competing order in a 
system, and the interactions between the two orders are 
strong, and a remarkably large and varied set of phase di- 
agrams are possible |?], |?], . The avoided critical point 
we have found here is a particularly striking example. In 
the context of stripe order, the situation is likely to be 
even more complex, since several distinct stripe orders 
have been observed in materials with stripy tendencies. 
Thus, when stripe order competes with superconductiv- 
ity, the full phase diagram should have multiple, possi- 
bly nested versions of the relatively simple phase diagram 
discussed in the present paper. 

There are a few additional observations we would 
like to make. The existence of vortex halos has many 
consequences [p9| , which we have not explored, for the 
character of the thermal vortex states. In particular, it 
gives rise to vortices that effectively have two||cJ core 
radii, £^ and £q, rather than a single vortex core radius 
as in conventional BCS superconductors. In this con- 
text, it is interesting to note that Wang et al^tt^, from an 
analysis of the Nernst effect in a variety of high tempera- 
ture superconductors, have recently adduced evidence for 
the existence of well defined vortex excitations at tem- 
peratures well above the superconducting T c . Moreover, 
from an analysis of the magnetic field dependence of the 



Nernst coefficient in La2-xSr ;r Cu04 and YBa2Cu307_5, 
they showed that there were two characteristic magnetic 
field strengths in the putative vortex liquid state, which 
they identified as corresponding to two distinct vortex 
core radii. If this interpretation is accepted, the smaller 
characteristic field, which they call B* , is to be associated 
with the maximal core radius, R* = W 4>q/2ttB* , where 
4>q = hc/2e is the superconducting flux quantum. At 
temperatures well below the zero- field T c , B* is typically 
in the range of 30T, corresponding to R* = 32A Since 
this length scale is comparable to the vortex halo ra- 
dius observed [Q in STM studies of near-optimally doped 
Bi2Sr2CaCu208+a, it is very tempting to associate R* 
with ^0. In the Nernst effect experiments, B* is seen to 
vary roughly linearly with temperature at low tempera- 
tures, so this identification could be tested by looking for 
similar temperature dependence of the vortex halo radius 
in STM experiments. 

Finally, we note that the picture advocated here has 
also a "dual" version Jl2[ |2^, in which one or the other 
form of stripe order is dominant, and superconducting 
order is sub-dominant. In this finite density of 

topological defects of the stripe state, dislocations, which 
could be induced by some form of shear or by disorder, 
plays the role of the vortex density the problem consid- 
ered in the present paper. Thus, a similar phase diagram 
as that in Fig. [j] (with the labels changed) can be con- 
structed for dislocation induced superconductivity in a 
stripe ordered phase. 
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